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Based on the model combining the spin double exchange and the lattice polaron, we have studied 
the colossal magnetoresistance phenomena observed in perovskite manganites Ri-^A^MnOs. First, 
effects of both the double exchange and the electron-phonon interaction on the transport property are 
investigated. We have evaluated the temperature dependent resistance and the magnetoresistance 
using the Kubo formula, and examined the crossover from tunneling to hopping regime of small 
polarons. Second, effects of the double exchange interaction on the lattice degree of freedom are 
^ — , explored. It is found that both the hardening of the phonon frequency and the reduction of the 

• phonon damping take place with decreasing the temperature. 
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I. INTRODUCTION 



The "colossal" magnetoresistance (CMR) manganites Ri-^A^MnOa (R= La, Pr, Nd; A= Ca, Ba, Sr, Pb) have 
recently attracted considerable attention due to scientific interest and potential applicability of their very large mag- 
netoresistance (MR) for 0.2 < x < 0.5 The most essential feature of their magnetic and transport behaviors is 
the existence of metallic conductivity and ferromagnetism. The magnetic transition at T c is closely connected with 
the resistivity peak at Tp corresponding to an insulator-metal transition (T c ~ Tp). The correlation between ferro- 
magnetism and metallic conductivity in Ri^A^MnOa was explained by Zener [Q in terms of the double exchange 
mechanism. There are mixed valent Mn ions (Mn 3+ and Mn 4+ ) as a consequence of hole doping by substituting 
R 3+ with A 2+ . In the double exchange model, conduction electrons in the partially filled e g levels of the d band are 
strongly coupled with the tightly bound d electrons in the t2 g levels by the on-site Hund's coupling, and mediate the 
ferromagnetic exchange interaction between the nearest neighbor S = | local spins formed from three d electrons in 
the core-like t 2g levels HH . 

Transport properties have been studied within the double exchange mechanism in favor of a magnetic polaron |l]-^j . 
Recently, Millis et al. || reported that the effective carrier-spin interaction involved in the ordinary double exchange 
Hamiltonian is too weak to produce the magnetic polaron effects. Instead, they suggested lattice polaron effects 
due to a strong electron-phonon interaction as a necessary additional extension [Q. They investigated a model of 
electrons, Jahn- Teller coupled to localized classical oscillators, within the dynamical mean field theory. Roder et al. 
jnj also examined the combined influence of the electron-phonon interaction and the double exchange on T c using 
the variational wave function techniques. But they could not treat the polaron transport and the lattice dynamics 
on an equal footing. In fact, the contribution of the lattice polaron to carrier mobility was pointed out earlier by 
Goodenough fl^| . 

There are many experimental evidences suggesting importance of the electron-lattice coupling in manganese oxides 
jljl 0. Near T c , dramatic changes are observed in the lattice degree of freedom - the anomalous lattice expansion 
beyond Griineisen law [O], and the shift of phonon frequency p5[-p7[, which all reflect that the lattice is closely 
related to the electronic and magnetic properties. However, detailed understanding of the interplay between the 
lattice dynamics and the electronic and magnetic properties remains to be resolved. 

In this paper, we have addressed two questions; i) what is the role of the electron-phonon interaction in CMR systems 
which arc known to have the double exchange interaction, and reversely, ii) how the double exchange interaction affects 
the lattice dynamics through the electron-phonon interaction. For these purposes, we first investigate effects of both 
the double exchange and the electron-phonon interaction on transport and magnetic properties. Employing the Kubo 
formula, we have determined the temperature dependent resistance and the magnetoresistance, and examined the 
crossover from a metallic tunneling state to an insulating hopping state in the small polaron transport. Second, 
to characterize the lattice dynamics in CMR compounds, we have considered the phonon degree of freedom in the 
presence of the double exchange interaction. We have studied the renormalization of the phonon frequency and the 
phonon damping constant. 

This paper is organized as follows. In section II, we present the model of conduction electrons coupled to phonons 
as well as the localized ionic spins in terms of the double exchange, whereby the temperature dependent resistance 
and the magnetoresistance are evaluated from the Kubo formula. In section III, we examine the double exchange 
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effects on the lattice degree of freedom. Finally, conclusions follow in section IV. Detailed calculational steps are given 
in Appendix. 



II. POLARON TRANSPORT 



Since Zener H] has proposed an interaction between spins of magnetic ions named " double exchange" , Anderson 
and Hasegawa B studied this mechanism in a system of Mn ions and a mobile electron with the transfer t between 
two Mn ions and the strong intra-atomic exchange integral J. When J is much larger than t, motion of the mobile 
electrons in Ri_ x A x Mn03 is described by the following double exchange Hamiltonian, 

ij 

where the hopping t%j connects neighboring sites, and fty is the angle between the directions of ionic spins at sites i 
and j. An exact quantum mechanical calculation gives cos^j- = ^gx^ , where S is the spin of a Mn ion and Sq is 
the total spin of Si, Sj and the conduction electron spin. In this study, we treat the double exchange part within the 
mean field theory following Kubo and Ohata jq], in which cos^j- is replaced by its thermodynamic average (cos^j-) 
determined by minimizing the free energy of the spin system. Then the propagation of an electron can be described 
as if it were moving in a mean field of highly disordered configurations of ionic spins. This approximation is known 
to work well at finite temperature, except for the very low temperature region (T ~ OK) where the spin dynamics 
becomes important. Within the present mean field theory, the double exchange plays a role, through ( S 2s+{ 2 ) = tC^Oj 
of increasing the bandwidth as T decreases below T c , accompanied by the ferromagnetic ordering (see Fig. 0(a)). 

In addition to the double exchange, the conduction electrons are scattered by the Mn-0 ionic motions in the MnC>6 
octahedra, which gives rise to a very strong electron-phonon interaction. The effective Hamiltonian incorporating the 
electron-phonon interaction is written as, 



H = Hcos-} ^ 4 +s c z + <"g«}a* + J2 ch^' R ' M q ( a ? + (2) 

iS q iq 

Here we adopt a model in which the single e g orbital is coupled to phonons assuming the electronically active e g 
band to be split, as in Roder et al.'s Jrl| . The present assumption is expected to be more effective if the model were 
generalized to include another physics such as the on-site Coulomb interaction, which might remove possible mid-gap 
states away from the Fermi level jlO| . 

The dc conductivity a can be obtained from the optical conductivity <j(u>) by taking the to —* limit, and <j(uS) can 
be determined by using the Kubo formula of the current-current correlation function, 

a(u) = — / dre-(4( T )j a (0)). (3) 
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Since the current operator J in narrow band systems is given by 
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a explicitly involves the four-site correlation function, 



J = it(co8-)e^2dct +s Cj, (4) 



^f^coAvEE^') l X dr(4(r)c j+s (r) 

S8> JJ> 



: ]'+S' c j' 



(5) 



In the isotropic case, the resistivity p corresponds to the inverse of a, p = 1/a. 

To evaluate a, let's consider the well-known polaron canonical transformation [fill; H = e s He~ s with S = 

— J2jq c j c j e ^ ' Rj Tf~( a q ~~ a -ij)- The transformed Hamiltonian is given by 

U = t(cos-) Y ct^CjX^Xj + J2 w g oJo? - A J2 c J c ii ( 6 ) 
jS q 3 
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with Xj = exp J^q e % ^' R i ^f-{a q — an( i A = X)<j 7T"- Inserting e s e 5 = 1 between each electron operator in 



Eq.(5), and using e b Cje 6 = CjXj and e b c^e 6 = c]Xj, one gets 

The intricate correlation function of Eq.(7) should be evaluated under the transformed Hamiltonian TL. Calculation 
can be further simplified with an approximation replacing the first term of Eq.(6) by f(cos|) z^2js(-^-j+S-^3') c j+5 c 3- 
This approximation is quite reasonable in view of that the coherent band-like motion of the polarons arises from the 
quantum mechanical tunneling between sites without changing the phonon numbers, which is governed by the matrix 
elements of (n< l \X'\ + gX j \riq) @. With increasing T, the polaron bandwidth decreases exponentially due to the term 

(xj +5 x,-)(=r(T)), 

r(T)=exp[-^Kf (iV 9 + l/2)], (8) 

where ug = (M q /uj q )(e l ^' s — 1) (see Fig. 0(b)). Under the above approximation, the mean field scheme of Kubo and 
Ohata M can begeneralized to include the phonon contributions which drastically reduce the magnetic transition 
temperature T c [pT[ , whereas the temperature dependent behavior of (cos|)(= j(T)) does not appreciably change. 

With the approximate Hamiltonian TL incorporating both "f(T) and T(T), the complicated four-site correlation func- 
tion in Eq.(7) can be disentangled into (Cj(T)cj+s(T)c^ l+s ,Cji) (xj(T)Xj + s(T)XJ l+s ,Xji) . These correlation functions 
and a can be evaluated in the straightforward fashion. Detailed calculational procedures are provided in Appendix. 
From Eq.(A12), the dc conductivity a is given as follows, 

ss ' 33' fcifca 

with a renormalized polaron band, tz = tq{T)T{T)^ s e~ lk ' 5 — A. Explicit expressions of ((j,j',5,S';T), £(T), and 
T)(j,j', 5, 5'; T) are given in Appendix. Keeping in mind that the auto-correlation function is most dominant for j = j' 
and 5 = 5', a is more simply obtained in the following form 

a = dt 2 7 (T) 2 e 2 Nz £ n %i {l - n £ Je ( ^-^ )/2T e - ( ^-^ )2/4C(T) e -«( T ) e ''( T )(VC(T)) 1 / 2 . (10) 

k±k 2 

Here z is the number of the nearest neighbors, and N q are the fermion and boson distribution function, respectively, 
and £(T), £(T), and r)(T) are also given by 

£(T) = ][> f | 2 (l + 2Ag, (11) 
C(T)=5»f| a [N^N. + l)] 1 / 2 , (12) 



V (T)^2j2\^\ 2 [N q (N q + l)} 1 / 2 . (13) 

Q 

Now let's investigate the qualitative behavior of p. One can carry out the numerical calculation of Eq.(10), assuming 
the simple square density of states (DOS) T>(e), 

. . N(l — x) W W . „, 

»(e) = — ~^-<e<—. (14) 
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Here the bandwidth W is given by wy(T)T(T) (w is the bare electron bandwidth and estimated to be w = 12|t| ~ 
2eV from the band structure calculation), and the Fermi energy tp is wj(T)T(T)(l — x) with x being the doping 
concentration. These treatments of DOS are based on the assumption that only the lower one of the split e g bands 
is active. We also assume that the most relevant phonon mode is the optical mode (~ ojq) which might be involved 
in the Jahn- Teller coupling. 

Numerical results for the resistivity are provided in Fig. ||. Two dashed lines represent resistivities of the polaron 
only model with given bandwidths of t up (= wT(T)) and t\ ow (— 0.75wT(T)). Here t up and ti ow correspond, respectively, 
to upper and lower limits of the double exchange factor 7(T). In both cases, p's exhibit peaks as a function of T, and 
the larger band width t up yields a smaller p with a peak at higher temperature. The resistivity peak as a function 
of T corresponds to the crossover from a quantum tunneling of the metallic phase to a self-trapped small polaron 
hopping of the insulating phase. Such features are characteristics of polaron systems, which are indeed observed in 
many oxide systems. In the high T limit, the resistivity has a thermal activation form of exp(A ff /T), characteristic 
of the semiconducting phase uM. Now taking the double exchange into account, the polaron bandwidth increases 
with decreasing T due to y(T^aJid accordingly the resistivity is given by the solid line (in Fig. ^J) with a peak at 
T c connecting the two polaron resistivity curves. This figure clearly demonstrates that the semiconducting behavior 
above T c is attributed to self-trapped lattice small polarons, and that the rapid fall-off in the resistivity below T c is 
attributed to the double exchange mechanism in addition to the lattice polaron effect. Thus, the combined model of 
the double exchange and the polaron provides a good description of the resistance anomaly observed in the experiment. 
The coincidence of the resistivity peak position Tp with T c originates from the mean field treatment of j(T) which 
neglects the fluctuation in the hopping of conduction electrons. 

Effects of the external magnetic field can also be taken into account in j(T) through the modified free energy 
due to the magnetic field. The behaviors of the MR's are shown in the inset of Fig. |^. With increasing the field 
intensity, the resistivity decreases and the peak position shifts to a higher T, and so the negative MR results. These 
results are quite consistent with the experimental observations. The MR peak Tmr is located near the resistivity 
peak Tp. In fact, Tmr,, Tp, and T c are the same in the present mean field treatment. It should be noticed that 
the magnitudes of the MR's in the figure are not large enough to explain the experimental CMR data quantitatively, 
suggesting that additional treatments might be required. One possibility is to incorporate the half-metallic nature 
of the ferromagnetic manganites |^,^lj , which is expected to suppress largely the spin-disorder scattering under the 
external magnetic field. 



III. LATTICE DYNAMICS 

As mentioned before, the phonon frequency becomes hardened as T decreases below T c jl5]-|l7|]. Interestingly, the 
phonon hardenings are observed for both optical and acoustic phonons in these systems. These frequency shifts are 
considered to be due to the change in the electron screening as T is lowered below T c . The hardening occurs in the 
metallic region (T < T c ), i.e., in the band- type tunneling regime where Rayleigh-Schrodinger perturbation theory is 
valid [Q. Therefore it is expected that the change of the bandwidth due to the double exchange factor "f(T) modifies 
the electron screening below T c . Note that the previous approximation for the canonically transformed Hamiltonian H 
corresponds to neglecting the phonon frequency renormalization, which seems to be too small to cause any appreciable 
change in the transport properties. In the metallic regime, the screening of the conduction electrons will be more 
easily described by the original Hamiltonian H of Eq.(2) rather than the polaron Hamiltonian H of Eq.(6). 

The renormalizcd phonon frequency uj q and the damping constants a q can be obtained from the following equation 

n 

{ui q -ia q j 2 =Lu 2 q - 2u q \M q \ 2 T{q, Cj q + iO + ), <I> g > a q , (15) 
where the electron screening function !F(q, u> q ) is given by 

The real part of both sides of Eq.(15) yields 

i fi ~* — 77,-* 

c .\ - 2. q \M q \ 2 -±- £ JL-t* (17) 

) g l k+q l k 
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It is important to note that the term J2 q ( n k ~ n fe+g) 7 '(%+<f — ' las ver y wea k temperature dependences (~ C + 
0(T/ Ep) 2 , -Ef being the Fermi level determined from tz). Hence the T-dependence of ui q comes dominantly from 
j(T), and Eq.(17) can be written as Cu q = /3/7(T)) 1 / 2 , where (3 contains all the T-independent terms. An explicit 

calculation of (3 is not available, but the order of its magnitude should be very small, 0(\M q \ 2 / (u) q Ep)) ~ 0(1CP 2 ). 
In Fig. ||(a), T-dependent behaviors of u)„ are plotted with respect to the external magnetic field strength H, and 
compared with the available experiment J16[ in the inset. It is seen that the frequency hardenings with decreasing T 
and with increasing H are qualitatively well explained. However, it is also apparent that some deviations exist between 
calculational and experimental results, particularly, near T c . These discrepancies might be ascribed to the mean field 
treatment of 7(T). Including spin correlation effects in calculating j(T) is expected to improve the agreement. It 
should also be noted that the phonon frequency hardening of Eq.(17) would be valid for both acoustic and optical 
phonon modes consistently with the experiments, because we have assumed the general form of the electron- phonon 
interaction in Eq.(2). 

Taking the imaginary part for both sides of Eq.(15), one gets the phonon damping parameter a q , 



2a a u) a = 2uj„\Ma\ 2 



1 _^ nj: — - 

m Im Z^ r; *._, 7j ,/JU-in+ - (18) 



99 i{T)""^ t %+i! - H -cj q H{T)-m- 

The imaginary part of the screening function is easily calculated by considering a parabolic electron band t% 

e £J~ \ t \s 2 k 2 , 



2a q = 2u q \M q \ 2 *- V{EF) ; ,19, 
2 v F q "f(T) 2 



where T>{Ep) and vf are determined from the parabolic band, |i|(5 2 /c 2 . In Fig. ||(b), T-dependent behaviors of a q 
are presented. Our results predict that the phonon damping parameter decreases with decreasing T, implying that 
the phonon is more sharply defined below T c . This feature in CMR systems is quite different from conventional 
observations of increased phonon damping parameter below T c for magnetic or strongly correlated systems. To our 
knowledge, no experimental reports are available yet on the phonon damping parameter. We think that the sound 
attenuation experiment will provide a better understanding of the nature of the electron-phonon interaction in CMR 
systems. 



IV. CONCLUSIONS 



We have extended the double exchange model to incorporate the strong electron-phonon interaction, and investi- 
gated transport and magnetic properties of CMR manganese oxides Ri_ x A x Mn03. We have found that the semicon- 
ducting behavior in manganites above T c is attributed to the effect of self-trapped lattice small polarons, and that 
the rapid fall-off in the resistivity below T c is attributed to the combined effect of coherent lattice polarons and the 
increased bandwidth via the double exchange mechanism accompanied by the ferromagnetic ordering of magnetic 
ions. Further, we have explored effects of the double exchange on the phonon degrees of freedom. The temperature 
dependent hardening of the phonon mode frequency observed in experiments is well described, and the reduction of 
the phonon damping constant is predicted below T c . 
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APPENDIX: CALCULATION OF dc CONDUCTIVITY a 



To evaluate the dc conductivity a, one should evaluate correlation functions of electrons and phonons under the 
Hamiltonian TL. Defining simply A and A' as 
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A = {e 1 ^ - 1)^, A' = e^ A 'i {e 1 ^' - 1)^, (Al) 

LO q UJ q 

X]{ T )X j+s {T)X] l+s ,X r is given by 

X}{T)X 3+S (T)X}, +s ,X r = n e -i(|A| 2 + |AT) e A*4e^ V A Qfe — (A2) 
tt • — Art -p - ^ 1 ? 7 " — A'*at — A'*at _ Art^p - ^ 177 " AA'%" 1W ? T 

Using e « e ? = e »e "« e , we see 

x}{r)X j+5 {r)x} l+5l X r — JJ^ e ~ T Hl A l 2 +l A 'l 2 ) e AA '* e_i " ,T e A * a f e ~ Aa ? ; (A3) 

<? 

where A = t±e~ luJqT — A'. Under the noninteracting phonon Hamiltonian, the thermodynamic average of (e X a ?e~ Aa «) 
is given by jl8j 

(e X * a te-^) = e~^, N q = (A4) 

Then (Xj(T)Xj + s(T)Xj, +s ,Xji) is obtained as follows 

{x]{r)X j+5 {r)x], +5 ,X f ) = exp[-$(i? J - 5', r)], (A5) 

$(i?, - R r ,6, 6', t)=J2 l%| 2 (! + 2iV «) 

- 2 £ ^0", ^)[JV 9 (JV 9 + l)] 1/2 co S K(r + (A6) 

9 

where U<r = {M q /uj q )(e 1 ^ - 1) and v^jj) = {M q /u) q )e^ R i {e 1 ^ - 1). 

The electron four-site correlation function (cj(r)cj+5(T)cj, +l5 ,c J v) is easily evaluated from the Hamiltonian H which 
transformed into the fc-space, 

(c](r)c j+5 (r)cl, +slCj ,) = £ » Sl (l - n^Je^^V^-^J'^-^Je*-^, (A7) 

where the renormalized polaron band is given as f7(T)(Xj +5 Xj) ^ 5 e lk ' s . 

The time integral of the correlation functions, Eq.(7) can be performed by the saddle point approximation. In the 
vicinity of the saddle point, the integrand becomes just a Gaussian, 

$(Rj - R r , 5, 6', r) ~ £(T) - ^(j, /, 5, 5'; T) + ((j, f, 6, 6'; T)z 2 , z = t + i/3/2, (A8) 

where £(T), i](j,j',S,S';T), and ((j,j',5,S';T) are, respectively, given by 

£(T)=£> 9 f (1 + 2^), (A9) 

9 

r,(j,j\S,6';T) = 2^ Uf (j, ^(j', S')[N q (N q + I)] 1 ' 2 , (A10) 

9 

C(j,j',S,S';T) = ^ 2 qVq {j,6)vt(j',5')[N q (N q + l)] 1 / 2 . (All) 

9 

^From Eqs.(7), (A5), and (A7), the dc conductivity a is evaluated and obtained as 
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"J 5 ' jj' fcife 2 

x e (*V*V/ 2 V (t V*V 2 / 4 W^ £ j». T)] i/2_ (A12) 
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FIG. 1. Bandwidths (hopping parameters in unit of t) as a function of the temperature, (a) Effect of the double exchange in- 
teraction, t^(T), (7(T) = (cosf }). (b) Combined effect of the double exchange and the electron-phonon interaction, tj(T)T(T). 
^2^\Uq\ 2 — 4.5 is taken. 



FIG. 2. The resistivity of the polaron only model (two dashed lines) and the combined model of the polaron and double 
exchange (solid line). In the polaron only model, t up and fi ow are the bandwidths corresponding to the the upper and 
lower limiting value of (cos|), respectively. Calculations are performed for the parameters, x = 0.3, w = 64T C ~ l.leV, 
uj a = 5T C ~ 0.08eV assuming T c ~ 200K, and ^-|w,j-| 2 = 6. The inset presents the resistance behaviors with H = 0T, 4.8T, 
and 9.6T. 

FIG. 3. (a) The phonon frequency shifts (u) q (T) — ui q ) /uj q for various magnetic field strengths, where uj q = ui q (l.lT c , H = 0T). 
In the inset, the shifts are compared with the experiment (under H — IT) for Lao. 7Cao.3Mn03 (Jeong et aZ.[16]). In the 
fitting, T c — 238K is taken from experiments. We have taken (3 — 0.07. (b) The phonon damping constants a q (T)/a q with 
a\ = a q (l.lT c ,H = 0T), are given. 
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